Cosmography of f{R) - brane cosmology 
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Cosmography is a useful tool to constrain cosmological models, in particular dark energy models. 
In the case of modified theories of gravity, where the equations of motion are generally quite com- 
plicated, cosmography can contribute to select realistic models without imposing arbitrary choices 
a priori. Indeed, its reliability is based on the assumptions that the universe is homogeneous and 
isotropic on large scale and luminosity distance can be "tracked" by the derivative series of the scale 
factor a{t). We apply this approach to induced gravity brane- world models where an /(i?)-term 
is present in the brane effective action. The virtue of the model is to self-accelerate the normal 
and healthy DGP branch once the /(_R)-term deviates from the Hilbert-Einstein action. We show 
that the model, coming from a fundamental theory, is consistent with the ACDM scenario at low 
redshift. We finally estimate the cosmographic parameters fitting the Union2 Type la Supernovae 
(SNela) dataset and the distance priors from Baryon Acoustic Oscillations (BAO) and then provide 
constraints on the present day values of f{R) and its second and third derivatives. 

PACS numbers: 98.80.-k,98.80.Es,11.10.-z 



I. INTRODUCTION 

The late-time acceleration of the Universe has been 
confirmed by several observations ranging from type la 
Supernovae (SNela) [l|, which brought the first evidence, 
to the cosmic microwave background (CMB) and the 
baryon acoustic oscillations (BAO) 0. More recently, 
gamma ray bursts (GRB), also if not properly standard 
candles, have been as well very useful at this regard 
They could, in principle, be useful to probe high redshifts 
with the aim to remove degeneracy of cosmological mod- 
els with respect to ACDM 0, 0] . While the recent speed 
up of the universe is a fact, we have yet no answer to the 
question: What is the "hand that rocks the cradle" ? 

If we assume that general relativity is valid on all the 
scales, even though it has been corroborated at most on 
the solar system range, then we require a component 
on the budget of the universe, that violates at least the 
strong energy condition to describe the current accelera- 
tion of the universe The simplest option at this re- 
gard corresponds to a cosmological constant, giving raise 
to the ACDM model which matches pretty well the obser- 
vations, but then we face the cosmological constant prob- 
lem. An alternative approach is to invoke a gravitational 
theory that deviates from general relativity on the appro- 
priate scales and at the same time being able to reproduce 
the big achievements of general relativity (cf. Refs. 
[Tsjl. The latter approach can be tackled in the context of 



brane- world models jl4l | , which are inspired in string the- 



ory, where our universe corresponds to a 4-dimensional 
hypcrsurfacc embedded on the higher dimensional space- 
time, usually dubbed the bulk. Several approach have 
been undertaken, for example in the context of induced 
gravity brane-world [l^, the self-accelerating brane 
of the Dvali-Gabadadze-Porrati (DGP) model is proba- 
bly the most famous [isj . 

The DGP model has gathered a lot of attention on 
the last years. As an induced gravity brane-world model, 
it contains two possible solutions, the self-accclcrating 
branch, which is asymptotically de Sitter, and the nor- 
mal branch. Despite this fact, the self- accelerating brane 
does not require any type of dark energy to describe a 
late-time inflationary period of the brane, it suffers from 
some theoretical problems like the ghost problem [l^; i.e. 
a degree of freedom that shows up when the brane is per- 
turbed and behaves on the brane effectively as a scalar 
field with the wrong kinetic energy. On the other hand, 
the normal branch is "healthy" in the sense that it does 
not suffer from the ghost problem but it requires some 
sort of dark energy to describe the late-time acceleration 
of the universe. 

In a previous paper [l9j . one of us proposed a mecha- 
nism to self-accelerate the normal DGP branch. More 
precisely, a generalized induced gravity brane-world 
model is proposed where the brane action contains an 
arbitrary f{R) term, R being the scalar curvature of the 
brane^. It is shown that an /(i?) (7^ R) term on the dy- 
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^ See Ref. |20|| for a brane-world model with an f{R) term in the 
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namics of a homogeneous and isotropic brane induces a 
shift on the energy density of the brane. This new shift 
term, which is absent in the DGP model, plays a cru- 
cial role to self-accelerate the generalized normal DGP 
branch of the model. In other terms, the generalized 
normal branch is asymptotically de Sitter without con- 
sidering any dark energy on the brane. 

In the present paper, we discuss the possibility to con- 
strain this model using a cosmographic approach [2l| . 
Cosmography relies on two crucial things: i) extract- 
ing the maximum amount of information from mea- 
sured distances, like the luminosity distances of SNela, 
ii) assuming that the universe can be modelled by a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) model 
on large scale without assuming a priori any dynami- 
cal theory to describe it. Now, why have we chosen this 
approach? for several reasons: i) for its simplicity. For 
example, the modified Einstein equation of the brane are 
of fourth order on the scale factor (due to the /(i?)-term 
in the brane action) and therefore very difficult to solve 
analytically. In the cosmographic approach we do not 
need to have an explicit solution for the evolution of the 
scale factor in terms of the cosmic time of the brane. ii) 
The approach is quite general in the sense that we do 
not have to specify which f{R) function we are dealing 
with. The only requirement is that /(i?) is an analytic 
function. 

The outline of the paper is as follows. In Sect. II, we re- 
view the model presented in jTof . In particular, we high- 
light how the model contains fixed points corresponding 
to de Sitter solutions (in absence of any matter on the 
brane); i.e. self-accelerating solutions. In Sect. Ill, we 
present the cosmographic approach we will follow. We 
write down all the quantities relevant of the model in 
terms of the cosmographic parameters. In Sect. IV, we 
constrain the model from a theoretical point of view, 
while Sect.V deals with observational constraints. Fi- 
nally, we summarize and discuss the obtained results in 
Sect.VI. 



the sake of simplicity, we have assumed a vanishing bulk 
cosmological constant, for a more general setup please 
see [l^- In addition, R is the scalar curvature of the 
induced metric on the brane, g, and k| is related to the 
Newtonian gravitational constant, G, through K4 = 8ttG. 
The function f{R) has mass square units. On the other 
hand, Cm corresponds to the standard matter Lagrangian 
of the brane. We recover the DGP model [l^, [lal when 
f{R) =R. 

From now on, we assume a homogeneous and isotropic 
brane with spatially flat sections. Therefore, the modified 
Friedmann equation can be written as 
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(2.2) 



The total energy density p is conserved and is given by 



p = Pm + Pf, 



(2.3) 



where 
Pf 



PmO 



^4 L 



3H\r--iRf-f) + 3HRr 



(2.4) 



where both energy densities pm and p / are conserved sep- 
arately. We will use the subscript to refer to quantities 
evaluated at the present time. The dot stands for deriva- 
tive with respect to the cosmic time of the brane and the 
prime for derivative respect to the scalar curvature of the 
brane. 

We are interested on the branch that generalize the 
standard DGP solution and therefore the modified Fried- 
mann equation (|2.2p reduces to 



H = 



6 



(2.5) 



II. AN /(i?)-TERM ON THE BRANE 

In this section, we review the model introduced in [l9} . 
The scenario corresponds to a 5-dimensional brane- world 
model whose action reads 



S = 




^5 ^'^4 

where is the 5D gravitational constant, is the 

scalar curvature in the bulk and K the extrinsic curva- 
ture of the brane in the higher dimensional bulk. For 



bulk. 



The other root of Eq. (|2.2p generalizes the Friedmann 
equation of the self-accelerating DGP solution. 

For latter convenience it is useful to rewrite Eq. (|2.5p 

as 



1 



-H = 



1^1 PmO 



Rf -f-QHRf") . (2.6) 



The parameter Tc = K'i^/{2K^ is the crossover scale. For 
1 ^ f'vcH, we obtain the Friedmann equation for 4- 
dimensional f{R) models. 

The Raychaudhuri equation for this model can be de- 
duced by taking the time derivative of Eq. (|2.6p , bearing 
in mind that the matter energy density is conserved, and 
it reads 



1 1 H 



1^1 PmO R^r + {R- HR)r 

2/' a3 2/' 



(2.7) 
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To obtain this equation we have as weh uscd^ R = 
6{2H^ + H). 

It can be shown that the brane contains fixed points 
corresponding to de Sitter solutions (once the matter con- 
tent is neghgible) [l^, therefore the brane enters a self- 
accelerating regime at some point along its expansion. 
In reference [l9[, it is shown what are the conditions to 
be fulfilled for the de Sitter solutions to be stable under 
homogeneous perturbations [l^. More precisely, we can 
associate an effective square mass to the perturbations 
and, as long as this quantity is positive, we can conclude 
that de Sitter solution is stable. 



III. COSMOGRAPHY 

A. General approach 

As we said, cosmography relies on the assumption that 
the universe is homogeneous and isotropic on large scale 
and no dynamical theory is assumed a priori [2l|. In 
particular it relies on the scale factor series expansion of 
a FLRW metric in terms of time [2l|; i.e. 

+b.Hiit-t,y^ + oiit-t„f) (3.1) 

where the standard cosmographic parameters are defined 

as m 



1 da 
H =- — 
a at 



1 d^a 



a change of sign of j (in an expanding universe) signals 
that the acceleration starts increasing or decreasing. 

Most importantly, the parameters {go, jo, 'o, so} can be 
used to evaluate different distances in the universe. This 
can be achieved by inverting the relation (j3.ip and bear- 
ing in mind that the distance, D, travelled by a given 
photon that was emitted at ti and detected at the cur- 
rent epoch is simply D = tg — ti (where we have set 
the speed of light to unity). Therefore, one can obtain a 
series expansion of the distance D in terms of the scale 
factor or redshisft, while the coefficients of the expansion 
are defined through the cosmographic parameters [2]| . 
The distance D can be related to several physical mag- 
nitude, for example the luminosity distance, the angular 
diameter distance and many more [2^ . These magni- 
tudes can be constrained observationally through SNela, 
BAO and, possibly, GRB data In fact, these data 
are useful to construct a cosmic ladder where any step is 
a cosmic indicator. Once the distances are constrained, 
we obtain as well constraints on the values acquired by 
the cosmographic parameter (see for example pll - [23l |). 
It is worthy to notice, at this regard, that given that the 
cosmographic approach is based on a Taylor expansion 
of the scale factor, or redshift, for data of GRB at high 
redshift (above z = 1), it is better to use the variable 
y = z/{l + z), introduced in [l^], instead of the redshift. 



B. Applying cosmography to f{R) brane-world 

In this subsection, we will relate the characteristic 
quantities defining the model introduced in Sect. II to the 
parameters {(?0i JO: 'oi so}- In addition, this will be done 
without specifying a particular f{R) model on the brane. 

We start reminding that the derivative of the Hubble 
parameter can be expressed in terms of the cosmographic 
parameters. Indeed, after some algebra, the following 
relation can be obtained: 



idh 



J = - 



dt^ 



H- 



(3.2) 



(3.3) 



s = -— H- 
a df^ 



H = H^{j + 3q + 2) 



(3.4) 



a dt^ 

These parameters are usually referred to as the Hubble, 
deceleration, jerk, snap and lerk parameters respectively 
(see [2l[ and references therein) . Their present day values 
(which we will denote with a subscript 0) can be used to 
characterize the evolutionary status of the Universe. For 
example, qq < denotes an accelerated expansion, while 



' We use Wald's book sign convention. 



H = H^ [s - Aj - 3q{q + 4) - 6] , 



(3.5) 



d^H/dt^ = H^[l-5s + 10{q + 2)j + 30{q + 2)q + 24] . 

(3.6) 

Now, the question is how our model can be 
characterized by these parameters, or, more pre- 
cisely, what can be said about the current values of 
f{Ro)J'iRo)J"{Ro)J"'{Ro)- In order to answer this 
question we have first to rewrite i?, i?, R, R in terms of 
q,j,s,l. This can be done with some algebra 
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i? = 6 + AHH^ , 



R = 6(H + AHH + AH' 



(3.7) 



Now using Eqs. ([SSI), and dSj]), we get 



i? = 6i/2(l-g) 



i? = 6H-^(j-g-2) 



6) 



(3.8) 



(3.9) 



(3.10) 



R = 6H^[l-s-2{q + 4)j-6{3q + 8)q-24] . (3.11) 

If we subs titute Eqs (EiSl), (EH), (ES]), (EH), (EH), 
([3^ . ([3?T0l) and ([3Tl|l in the Friedmann and Ray- 
chaudhuri equations, and evaluate them at the present 
time, we could obtain, in principle, the current val- 
ues of f{Ro),f'{Ro)J"{Ro)J'"{Ro)- However as we 
have only two equations, the Friedmann relation and 
the Raychaudhuri equation, we require more informa- 
tion to define completely the model. At this respect, no- 
tice that the effective gravitational constant on the brane 
GeS = G/f (see the Friedmann equation ()2.6p ). there- 
fore we can assume, as a prior, that /'(i?o) ~ 1 such 
that the current value of the gravitational constant coin- 
cides with the Newtonian one. Further information can 
be obtained through the equation satisfied by H- At this 
respect, we take the time derivative of Eq. (|2.7p and we 
obtain 



H 



1 [HH - H^)f' - HHf'R _ 

q7c WW ~ 

R^f" + {R- RH)f" + Klpmpg-^ 
2/'2(/"7?)-i 
+ {3RR - Hi{^)f"' 
V' 

(R-RH- RH)f" - SKlHp^oa-- 
2/' 



(3.12) 



where f^™^ = d^f /dR^. As third assumption, we take 
into account the power series 

/(i?) ^ f{R^) + f'{Rp){R-Ro) + \f"{Ro){R-Ra? 



+ -f'\R,){R-R,)\ 



i.e. at low redshift, the /(i?)-function is well approxi- 
mated by its Taylor expansion up to the third order"^. 

Now we can finally substitute Eqs. p.3p . p.4p . 
(|331) . dSH), daH), (|XTO| and (f3TI|) in the Fried- 
mann constraint p.6p . the Raychaudhuri relation ()2.7p 
and the complementary equation p.l2p . We evalu- 
ate them at the present time. Notice that Eqs. (|2.6p 
and (|2.7p can be expressed as linear combinations of 
f{Ro)J'{Ro)J"{Ro)J"'{Ro) at z = 0. This is not the 
case for equation p.l2p as it is quadratic on /"(i?o)- So, 
we will proceed as follows, we obtain f{Ro) as a linear 
combination of /"(i?o) using Eq. (|2.6|) . 



fiRo) = 6 



(r!,„ - 1) i/2 + ^ ( i?o - GHoRof ] - -Ho 



1 /„ ... ^ „„\ 1 
6 V r, 

(3.14) 

where Hq = K\pma / {iH'^) . Similarly, we can write 
f"'{Ro) as a linear combination of /"(i?o) using Eq. (|2.7p . 
i.e. 



f'iRo) 



iHl^ra + -ffo(2 + ^) + (i?o - HRo) f'iRo 



Rq 



(3.15) 



Then we rewrite Eq. p.l2p as follows 

a2f"iRof + ai f'iRo) + ao = 0, (3.16) 

where 

02 = RoiRo-RoHo), (3.17) 
fli = Rof"'iRo) + '^H^QmRa — (^iio — Rq — RqHo^ 

+ ^ (3.18) 

ao = — (sRoRo — HoRo^ fiRo) + ^H^D.^ — 2Ho 
Ho Ho — Hq 



rcHl 



(3.19) 



Even though the previous equation looks quadratic in 
f'iRo), it is not the case because ai is a linear function of 
f'iRo) and therefore this term contributes quadratically 
in /"(i?o)- Once we substitute Eq. ((XT51) on Eq. ((XTCl) . 
we obtain a linear equation for f'iRo). 

Finally, we obtain the following results: 



, (3.20) 
, (3.21) 
, (3.22) 



fiRo) 


Aon„,+Bo+CoircHo)-^ 




V 


f'iRo) 


A2n„,+B2+C2ircHo)-^ 


(6i72)-i 


V 


/"'(i?o) 


A3nrn+B:i+C:iir,Ho)-^ 


(6F2)-2 


ijo - go - 2)P 



(3.13) 

^ In what follows wc assume f'{Ro) = 1 and f^'^^{Ro) — 0. 
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where Ai,Bi,Ci and V with i 
q, j, s, I which are defined as 



0,2,3 arc functions of 



^0 



Bo 



Co 



(jo - go - 2) lo 

(3so + 7 Jo + GQo + 41go + 22) sq 

[(3(?o + 16) 30 + 20ql + 64qo + 12] jo 

{3q^ + 2bql + mql + 12qo + 20) , (3.23) 

- {joqo -ql- 2qo) Iq 

[3goSo + (4(?o + 6) jo + Gqo + 44(7o + 22go 
-12] so 

[2j2 + (3q2 + lOqo - 6) .70 + I7ql + 52q^ 
+54go + 36] jo 

3q^ + 28qf, + 118q^ + 72ql - 76qo - 64, (3.24) 

- {ja - qo-2) lo 

[3so + (3qo + 1) jo + 'iql+ 41qo + 34] so 
[il - (9o - 90 - 6) jo + 'oql + A'6ql 
+50(70 + 4] jo 

(q^ + 3g3-80q2-144go- 68), (3.25) 



A2 
B2 



9so + 6jo + 9ql + 66qo + 42, (3.26) 

- {6 {qo + 1) So + 2 [jo + {qo - 1)] Jo 

6q^ + 50q^ + 74qo + 32} , (3.27) 

- {3 (1 + go) so + [jo - (Qo +qo + 2)] jo 

4g3 + 29ql + 42go + 18} , (3.28) 



-43 = -3 [lo + so - 3 {qo + 4) jo - 15go 

-26go - 4] , 
B3 = 2[(l + go)^o + (<7o+jo)so 

- (jo + 2go + 6go + 3) jo 

- {I5q^ + A2ql + 39<zo + 12)] , 
C3 = (1 + qo) k + {jo - 90 - 90 - 1) So 

- (.70 + 9o + 4go + 2) Jo 

- (<z;; + 26g3 + 69gg + 64qo + 20), 

V = - (.70 - go - 2) / + (3so - 2jo + 6<Zo' 
+ [(3go + 10) J + llgg + 4go - 18] jo 
+ 3q^ + Mqf, + l&Oql + 246go + 104. 



50go 



(3.29) 



(3.30) 



(3.31) 
40) so 

(3.32) 



We have split the expressions of /(i?o), f"{Ro) and 
/"'(i?o) into three pieces involving the functions Ai^Lrm 
Bi and Ci{rcHo)~^ , where Ai, Bi and Ci are defined exclu- 
sively in terms of the cosmographic parameters. The first 
term Aiflm account for the contribution of matter to the 
/(i?)-function^. The second one Bi is a purely geomet- 



* It is worth noticing that wc are developing our considerations in 
the Jordan frame so the standard matter is minimally coupled 
to the geometry. 



rical one. The third one takes into account the effect of 
the extra dimension; i.e. it involves the crossover scale re- 
Not surprisingly, if we switch off this term: i.e. 1 ^ rc, we 
recover exactly the results obtained in [2l| corresponding 
to a standard 4-dimensional f{R) scenario. 

In summary, for a given set of values of the cosmo- 
graphic parameters we can deduce the function f{R) 
through the expression p.l3p . Notice that the opposite is 
not possible because the equations (3.23)-(3.32) are non- 
linear in {go, So, ^0, jo}. Moreover, by specifying a given 
function f{R), we do not obtain a unique evolution for 
the brane because the modified Raychaudhuri equation 
is of fourth order in the scale factor. 



IV. PARAMETERIZING THE 
COSMOGRAPHIC PARAMETERS 

In order to get a first hint on the possible values of 
f{R) and its derivatives we adopt the following strat- 
egy: the cosmographic parameters will be calculated for 
a given dark energy phenomenological parameterization. 
The best and simplest one is the ACDM model. Next, we 
will evaluate those parameters using the recent data of 
WMAP7 and the constraint on the crossover scale Tc (see 
[25} for details). Through these results, we can constrain 
the f{R) function as we will show below. This is a mini- 
mal approach but it is useful to probe the self-consistency 
of the model. 

The cosmographic parameters for the ACDM model 
read 



Ho 
H 



i i Lf-,i — 17- 

2 a-^ 



Ho 
H 

Ho 
H 



Ho 
H 



1 

1 - 2^,, 

2 a3 

1 - 2fi„ 



' ^ 2~^ 

7 n"^ 

2 



(4.1) 
(4.2) 

(4.3) 



35 



which, evaluated at the present time, give 21 1 

3 



10 

jo 

So 



-1 



1, 

9 

1 - ^Jl. 



27 



(4.4) 

(4.5) 
(4.6) 
(4.7) 

(4.8) 
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Inserting the previous equations in the equations 
(3.22)-(3.31), we obtain 



^0 

Bo 
Co 
A2 
B2 
C2 

-43 

V 



63 



27 



243 



(4.9) 



81 72Q 

§^f„ + f^L(4.io) 



63172„ 
63, 



81 



16 



81. 



441 , 27 , 

18917„-— 172„-84--f)^, 

4 

8 ™ 16 ™ 



(4.11) 
(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 
(4.18) 



63rJ^ ■ 4 --™ ' 16 

It can be checked that if 1 ^ rc] i.e. in absence of an 
extra dimension, the function f{R) reduces to f{R) ^ 
R-2A because f"{Ro) = and f"'{Ro) 0. This can 
be assumed as a consistency check. However, as soon as 
the effect of the extra dimension is switched on, i.e. Tc 
is finite, the coefficients Ci with i = 0, 2, 3 play a crucial 
in defining the shape of the function f{R). Indeed, we 
obtain 

Aon,n + Bo 



^GRo 

Fgr, 
Fgfi.3 

Fig, 



Ro - 2A 
= 0, 



= 0, 



Fja„ — — 



V 

+ 1 = 

A2^m + B2 

V 

{]0 -qo- 2)P 
CojrcHo)-^ 
V 

-112- 18a„ + 91}^ 

112 + mn,r, + 27nf^ 

C2ircHo)~^ 
V 

4 -25217„ + 147»,^„ + 112 + 18»f„ 
3 n^(112 + 6017„ + 27172„) 

CzircHo)-^ 



{rcHor\ 



Uo - 90 - 2)2? 

26 + 3Qr. 



112 + 60an + 27f72 



■{rcHoY 



4.19) 
4.20) 

4.21) 
4.22) 

4.23) 
4.24) 

4.25) 

r\ 

4.26) 
4.27) 

4.28) 



For clarity, we have split the right hand side (rhs) of 
Eq. (3.20) into two pieces: Fgrq and F/g^, the first one 



takes into account the pure relativistic contribution while 
the second one takes into account the effect of the extra 
dimension. A similar procedure has been followed with 
the rhs of Eqs. (3.21) and (3.22). 

We consider the following observational conservative 
values ftm = 0.266 and ilr^ = 10^"* where ilr^ ~ 
{4:rcHQ)~^ 0, [2^ and we obtain the values reported be- 
low: 



FgRo 
Fgr2 
FgRs 

FiGo 
F1G2 
F1G3 



0.867, 

0, 
0, 

-0.018, 
0.161 , 
0.003, 



with the errors evaluated as in [3, [11] . The previous re- 
sults show that, although FjGi are different from zero, 
they are relatively small in comparison with the present 
day main contribution Fgr„] i.e. the standard relativis- 
tic term. In summary, the model deviates just slightly 
from the pure ADGP model^ [Mill- This smafi devia- 
tion is enough to obtain self-acceleration without invok- 
ing any kind of dark energy contribution on the brane. 
On the other hand, if a similar analysis is carried out 
for a given f{R) function in a 4-dimensional model, it 
turns out that the /(i?)-term match completely that of a 



Hilbert-Einstein action plus a cosmological constant [21 



Most importantly, we see that the model we have ana- 
lyzed is consistent with the ACDM model because the 
cosmographic parameters of the ACDM can be matched 
to those of an f{R) branc-world scenario. 



V. OBSERVATIONAL CONSTRAINTS 

In order to constrain the model, i.e. to estimate the 
function f{R) through its own value and that of its 
derivatives at the present time, we need to constrain ob- 
servationally the cosmographic parameters by using ap- 
propriate distance indicators. Moreover, we must take 
care that the expansion of the distance related quantities 
in terms of (go, Jo, so,lo) closely follows the exact expres- 
sions over the range probed by the data used. Taking 
SNela and a fiducial ACDM model as a test case, one 
has to check that the approximated luminosity distance^ 
deviates from the ACDM one less than the measurement 
uncertainties up to z ~ 1.5 to avoid introducing any 
systematic bias. Since we are interested in constrain- 
ing {qoija, So, ^0), we will expand the luminosity distance 
Dl up to the fifth order in z which indeed allows us 



The ADGP model corresponds to the normal DGP branch en- 
dowed with a cosmological constant and filled with matter. 
® See [2]| for the analytical expression. 
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X XBF jx) Xmed 68% CL 95% CL 

h 0.744 0.750 0.750 (0.725, 0.775) (0.701, 0.802) 

go -0.43 -0.44 -0.45 (-0.48, -0.41) (-0.51, -0.36) 

jo -0.35 0.01 0.01 (-0.11, 0.14) (-0.33, 0.35) 

so -1.3 0.4 0.4 (-0.3, 1.0) (-1.2, 1.8) 

lo 14.7 -0.6 -1.0 (-4.6, 3.7) (-11.3, 11.7) 



TABLE I: Constraints on the cosmographic parameters by 
jointly fitting the Union2 SNela sample and the BAO data. 
Columns are as follows: 1. parameter id; 2. best fit; 3., 4. 
mean and median from the marginalized likelihood; 5., 6. 68 
and 95% confidence ranges. 



to track the ACDM expression with an error less than 
1% over the full redshift range. We have checked that 
this is the case also for the angular diameter distance 
Da = Dl{z)/{1 + zY and the Hubble parameter H{z) 
which, however, we expand only up to the fourth order 
to avoid introducing a further cosmographic parameter. 

In order to constrain the parameters {h, go, ia, sq.Iq), 
we use both the Union2 SNela dataset [2^ and the BAO 
data from the analysis of the SDSS seventh release poj . 
We then consider the following likelihood function : 



C{p) = CsNelaip) X CsAoip) 



(5.1) 



where p is the set of model parameters and we have de- 
fined the likelihood function for the probe i as : 



A(P) 



(2^)A^./2|C,|l/; 



■ exp 



(5.2) 



For SNela, AsNeia is MsNeia (with NsNeia = 557) col- 
umn vector with elements computed as : 

AsATe/aj = fJ-obs{Zj) - Hth{Zj,p) (5.3) 

fithiz) = 25 + 5\og Dl{z, p) , (5.4) 
while the CsNeia is a diagonal matrix. For BAO, we set : 



dthiz, p) 



^BAO 



^(Zd) 



J = dobs{z) - dth{zj,p) 



Dv{z,p) 



rs{zd) 



{l + zfDl{z,^)cz 
H{z,p) 



(5.5) 



-1/3 



(5.6) 



where we set the sound horizon distance to the drag red- 
shift as rs{zd) = 152.6 Mpc. Percival et al. [1^ provide 
estimates of dz for z ~ (0.20, 0.35) and the correspond- 
ing covariance matrix that we use as input in Eq. (j5.5p . 
We remember the reader that we use a fifth order ex- 
pansion in z for both Dl{z) and Da{z), while H{z) is 
expanded to the fourth order only. Since the BAO data 
are at low redshift, the resulting approximated expres- 
sion for dth{z) closely follows the exact values. Finally, 
we also use a Gaussian prior on h from local distance 
measurement so that (j5.2p reduces to a Gaussian centred 
on h = 0.742 and with variance ah = 0.036 [30| . 

In order to sample the five dimensional parameter 
space, we use a Markov Chain Monte Carlo algorithm 
running two chains (with 125000 point each) and check- 
ing the convergence according to the Gelman - Rubin cri- 
terium (R—l < 0.1). The resulting constraints are sum- 
marized in Table I where we give the best fit parame- 
ters and the constraints over the single pi obtained by 
marginalizing over the other ones. As a general remark, 
we find that these constraints are i n ag reement with pre- 
vious constraints in literature (23l . l3l| . Note, however, 
that our confidence ranges turn out to be narrower than 
usually found. This is likely due to our inclusion of the 
lerk parameter Iq. In a sense, we are now better ap- 
proximating the (unknown) actual distances and Hub- 
ble parameter so that not all the possible combinations 
of (/i, go, jo, So) are possible, but only the ones that are 
compatible with the constrained lo- 
in order to translate our constraints on the cosmo- 
graphic parameters on similar constraints on f{R) and 
its derivatives, we should just use Eas. p.22p - p.32p eval- 
uating them along the final coaddcd and thinned chain 
and then looking at the corresponding histograms. To 
this end, however, we should set also the values of 
and fi^c (and hence VcHq = l/2-y/iX~) which are not 
constrained by the fitting analysis described before. To 
partially overcome this difficulty, we adopt the following 
strategy. Defining for shortness 



/o = 



/GRo) 

6i72 



/2 = 



f"{Ro 



h = 



f"'{Ro 



we first constrain these quantities setting fi,.^ = lO^'' 
and varying Q.m along the chain using Q,m = uiMh~'^ 
with the physical matter density um = 0.1329 in agree- 
ment with the WMAP7 data. Note that we are neglect- 
ing the uncertainty on w^/ since it is much lower than 
those on the cosmographic parameters. We also stress 
that, although the fiducial value for lum has been ob- 
tained for a ACDM model, it should be unchanged for 
any model which reduces to the GR -|- matter domina- 
tion at the CMBR epoch as is our case. We can then 
scale the results to a different value of r^-ffo noting that, 
by simple algebra, we get from Eq. (|3.22p : 



-0.55 -0.50 -0.45 -0.40 -0.35 -0.30 
'/u 



-0.55 -0.50 -0.45 -0.40 -0.35 
'/u 



-0.55 -0.50 -0.45 -0.40 -0.35 -0.30 



FIG. 1: Isolikelihood (68, 95 and 99% CL) contours for the fit to the SNela and BAO data. In each panel, we marginahze the 
given cosmographic parameter with respect to qo. 



X 


xbf 


(x) 




68% CL 


95% CL 


fo 


0.897 


0.912 


0.912 


(0.876, 


0.949) 


(0.828, 0.992) 


h 


0.126 


0.163 


0.161 


(0.140, 


0.185) 


(0.116, 0.220) 


h 


-0.130 


-0.139 


-0.142 


(-0.181, 


-0.101) 


(-0.240, -0.004) 


ao 


-0.0190 


-0.0168 


-0.0167 


(-0.0180, 


-0.0155) 


(-0.0201, -0.0141) 


Po 


1.0190 


1.0168 


1.0167 


(1.0155, 


1.0180) 


(1.0141, 1.0201) 


OL2 


0.0130 


0.0190 


0.0191 


(0.0172, 


0.0209) 


(0.0137, 0.0236) 


P2 


0.9870 


0.9810 


0.9809 


(0.9791, 


0.9828) 


(0.9764, 0.9863) 


Q3 


3.1091 


0.0071 


0.0100 


(0.0037, 


0.0140) 


(-0.0272, 0.0263) 




-2.1090 


0.9929 


0.9899 


(0.9860, 


0.9962) 


(0.9736, 1.0272) 



(5.9) 



TABLE IL Constraints on the fiducial /; values and on the 
scaling coefiicients {ai,j3i) from the Markov Chain for the 
cosmographic parameters. Columns are as in Table L 



The constraints on the fiducial fi and the scaling pa- 
rameters (ai, Pi) obtained by evaluating these quantities 
along the Markov chain for the cosmographic parameters 
are summarized in Table IL Considering the median 
values and the quite narrow confidence ranges, we find 
that that the scaling parameters (ai,/3i) are well con- 
sistent with the /; being linear functions of the inverse 
of the crossover scale Tc hence allowing us to easily 
estimate the impact of uncertainties on this parameter 
on the final estimate of the present day values of f{R) 
and its derivatives. Somewhat surprisingly, the fiducial 
fi are reasonably well constrained notwithstanding the 
large uncertainties on the cosmographic parameters. 
Such a result can be qualitatively understood noting 
that fi depend on (<?o, Jo, sq, ^o) through a ratio of 
coefficients so that it is possible that a variation in the 
numerator is compensated by a similar variation in the 
denominator in such a way that the final fi is unaltered. 
As a consequence, the dependence on the cosmographic 
parameters is made weaker thus reducing the impact of 
the parameters uncertainties. 



VI. CONCLUSIONS 



fid 



ai ^ h ft = Qfi 



TcHq 



1/2 



+ ft (5.7) 



with the quantities labelled fid are obtained for the fidu- 
cial ri^c value and we have defined (for i = 0, 2, 3) : 



(5.8) 



Cosmography is a useful method to give a picture of 
the observed universe considering minimal assumptions 
(isotropy, homogeneity, Taylor series expansion of dis- 
tances) without choosing any dynamical model a priori. 

In this paper, we have taken into account the problem 
to test brane-cosmology, where an /(i?)-term is present in 
the boundary 4D-action, by cosmography. Being ACDM 
a realistic picture of the today observed universe, we have 
adopted ACDM observational results as priors for our ap- 
proach. Wc assumed the f{R) function to be analytical 
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in order to evaluate the higher-order curvature contri- 
butions with respect to general relativity contribution, 
i.e. f{R) — R. The results arc encouraging since small 
higher-order deviations with respect to general relativity 
give dynamical behaviors, consistent with observed cos- 
mic acceleration, without introducing dark energy terms. 

However, the approach should be consistently probed 
at small, medium and high redshift by selecting suitable 
standard candles or, at least, reliable distance indica- 
tors at any scale. Despite of this technical difficulty, 
the method outlined here deserves further investigations 
since it is connecting a fundamental theory, as the DGP- 
brane model, with data coming from precision cosmology. 
We have here addressed this point in a preliminary way 
by only using SNela and BAO, but other probes (such as 



GRBs) may be added to further narrow the constraints 
on the present day values of f{R) and its second and 
third derivatives with respect to R. 
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